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' In a molecular semiconductor, the carrier is dressed with a polarization cloud that 

we treat as a quantum field of Frenkel excitons coupled to it. The consequences 
of the existence of this electronic polaron on the dynamics of an extra charge in a 
, material like pentacene can thus be evaluated. 
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Molecular organic semiconductors in general and oligoacenes in particular are seriously 
being considered today as materials for application in optoelectronic devices. Pentacene, 



for example, has already been used with success in organic field effect transistors 0, 0, HI- 



Critical to their eventual use as electronic materials, therefore, is a proper understanding 
of the charge transport mechanisms in these molecular solids. Currently, the paradigm in 
force describing the electronic properties of oligoacenes views carrier transport in terms of 
lattice polaron theories applied to narrow band crystals 0, El IE 00- However, we believe, 
as other authors [9|, [l0| have also pointed out, that this approach omits an important effect 
specific to organic materials composed of large polarizable molecules, i.e., adding a carrier 
onto a molecule in such a crystal, creates a polarization cloud around the carrier, also-called 
electronic polaron or Coulomb polaron [llj], which accompanies the itinerant charge as it 
moves about in the solid. This effect has actually been confirmed by photoemission experi- 



ments [Tj, [13( wherein the spectra of individual molecules in the gas phase are compared to 
the corresponding spectra in the solid phase [l2| or within a molecular cluster [l3[. 

In these crystals, the time scale for establishing molecular polarization is related to the 
energetic distance e between the ground state and the first excited state of the neutral 
molecule, generally on the order one electronvolt. The characteristic time, given by t £ = h/e 
(where h is the Planck constant), is thus on the order of 2.10 -15 s. This is a very fast 
process, about ten times faster than the characteristic time for carrier motion from molecule 
to molecule tj = h/J, where J ~ 0, 1 eV is the bare transfer integral between two adjacent 
molecular sites along the easiest direction of propagation in the crystal [8]. On this basis, 
one would thus expect the Coulomb polaron to readily follow the extra charge in its motion 
through the crystal without any significant scattering or friction. The general feeling among 
physicists therefore is that polarization effects are minor and that electron-phonon processes 
probably explain the bulk of transport phenomena in these materials. 

We will show, however, that the electronic polaron does actually play an important role in 
the dynamics of an extra carrier, in particular, by (i) significantly renormalizing the transfer 
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integrals in the perfect crystal, (ii) enhancing the effect of any type of disorder including 
thermal disorder, and (Hi) creating correlations between diagonal (on-site) disorder and non- 
diagonal (transfer integral) disorder. In contrast to early calculations (reviewed by Silinsh 
and Capek (HI), as well as more recent ones fTo| . which use a self-consistent approach in the 
limit of vanishing molecular overlap to compute the polarization energies, we will present 
here a theoretical approach which accounts for the finite bandwidth effects by treating the 
entire solid quantum-mechanically. This is essential for studying the transfer of the itinerant 
carrier through the solid. The treatment of polarization here is applicable to all types of 
organic solids composed of large conjugated molecules, disordered or not and provides good 
insight into the possible roles of polarization in charge transport. 

To demonstrate our points, we start from a model hamiltonian, inspired by the pioneering 
work of Capek on the polarization of a molecular trimer We will consider both the tight 
binding charge propagation from site n to the neighboring sites n + h and the consequent 
polarization of the molecular sites I coupled electrostatically to the extra charge. Thus, the 
annihilation and creation operators a and d + describe the narrow band J excitations while 
the operators b and b + take into account the internal electronic degrees of freedom in the 
neutral molecules responsible for polarization. 

When the molecule is anisotropic, like pentacene, the polarization excitations are direction 
dependent. It is convenient then to describe each molecule I located at position 17 in the 
frame of the principal axes of its polarizability tensor at. In this ortho normal frame where 
the unit vectors can be called (e),e|,ef), the dipolar transition tensor fii is also diagonal 
and can be represented by the three components //^ 

Hi = e(</>ffM0l) (1) 

where 4>h is the non-degenerate HOMO state and 0^ the particular LUMO state excited 
by a dipolar transition in the direction i. In the spirit of a tight binding model, the <p l L are 
three different linear combinations of the it orbitals selected when the particular direction % 
of the dipolar transition is considered. These states have essentially the same energy e with 
respect to the HOMO level [la ]. In this approach, the polarization of the molecules I results 
from a small depopulation of the HOMO in favor of the LUMO. This excitation process, 
internal to the molecules is represented by the operator pairs h\. and b~\.. 

The quantum operator representing the dipolar moment on each site I is the vector 

3 

d e = 2J de, el with d iz = ^ {b t + bj)i (2) 

i=l 

The diagonal polarization tensor ae is then written as (at)u = 2 fi^/e. 

With this notation, the hamiltonian for the extra charge, in the presence of the polariza- 
tion cloud can then be written in the dipolar approximation as: 
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where J n , n +h are the transfer integrals between sites n and n + h and 



(We, 
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are the dipolar interactions between pairs of molecules in the structure, responsible for 
the Van der Waals contribution to the cohesive energy. 

In fact this hamiltonian describes a simplified version of the motion of a charge carrier 
interacting with a quantum field of Frenkel excitons of energy e [11]. 

In molecular crystals like acenes, the motion of the excess charge is slow compared to the 
relaxation time necessary for the polarization of the electronic orbitals of the molecules sur- 
rounding the charge carrier, i.e. J <C e. To approximate ground state of the hamiltonian |31 
we have thus used a variational method with a dressed carrier trial function. A perturbation 
theory would also give such solution when r £ <C tj. The eigenstates of the hamiltonian are 
represented as a superposition of local states, taken to be the product of a local electronic 
state |n) and the polarization state \x{ n )) of the surroundings, associated with the carrier's 
occupation of site n : 

n 

The local polarization states which constitute \xi n )) can be expressed in terms of the 
unitary translation operators Ue where 1 0^) is the ground state of the £-th molecule and 

\x{n)) = ®l[Ue(n)\0t) (6) 

t 

U e (n) = exp (]T fain) K - X^n) &+' 
\i=i 

The unitary operator Ug(n) represents a translation of the molecular state due to the 
electric field of the charge placed on site n: 

Ur\n) h Ue(n) = K - X tfi (n) I t (7) 

U^{n)blUe{n) = 6+ - X^n) i t 
where Ie is the identity operator. 

The functions X^j(n) are determined by minimizing the expectation value of the Hamil- 
tonian El with the variational wave function |3] and El above : 

3 
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where we have used the Weyl identity 

(0,| e ( sh t- s+h *) e(^'-^ + ) \0 e ) = e (^-(l*l a+ W a )/ 2 ) (9) 

We now minimize the energy of equation El with respect to the variational functions, the 
X^j(n)'s. Since we are concerned with the stability of a slow carrier, we need only minimize 
the potential-like energy. Thus to the lowest order in J/e, we get 

e*£,<(n) + ^ ^Z~-rf + E E (*W<J (XmM + X* mJ (n)) = (10) 

' ' m j=l 

Substituing ITU1 in 19*1 the local polarization energy becomes: 

ifn i=l 

Within the subspace of states \ip) representing the dressed-electrons defined by relation 
El , the total energy minimum El can thus be written as: 

E = 22E p (n) \u n \ 2 - 2j J n ,n+hU* n+h U n (12) 
n n,n+h 

Then, the hamiltonian El describes a dressed carrier (or a quasi-particle) at site n having 
a potential energy E p (n) and an effective transfer integral 

T - T p-SoO) 

"n,n-\-h < - , n,n+/i c 

with 

s »( ft ) = \ E E - *y» + ; Ep(0) t £p(ft) (13) 

£ i=l 

where 7 is a coefficient of order unit. 

The effective quasi-particle hamiltonian H with parameters renormalized by the Frenkel 
exciton field is just the tight-binding hamiltonian : 

H = ^E p {n) \n)(n\ - J n , n+h \n)(n + h\ (14) 

n n,n+h 

where E p (n) and J n ,n+h are defined by relations ITT1 and IT^l 

In the case of a perfect crystal, the polarization energy E p is uniform and shifts the 
ground state everywhere by about 1 eV as observed in acenes jl^, EH- Correlatively the 
bare bandwidth is significantly narrowed (independently of the temperature). Values of E p 
(from and J/ J have been reported in table HI We can conclude therefore that at low 
temperatures in a perfect molecular crystal, a band model with an appropriate renormalized 
bandwidth describes the extended ground state properly |l4j (claim (i) above). 

The situation changes greatly with disorder, however small, enters the system. The local 
changes of polarization energy E p due to disorder are greater than any other electronic pa- 
rameter change, transfer integral, HOMO and LUMO positions, etc. In particular, thermal 
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disorder arising from librations or low energy intermolecular phonons is interpreted as being 
static on the time scale Tj and induces non negligible polarization energy variations. 

Indeed, consider any kind of static (or thermal) disorder which can be described by a 
random variable in the lattice. The polarization energy E p is in turn also random. We 
denote S n the fluctuation of E p from site to site. Then J n ,n+/i a l so becomes random through 
the fluctuations of Xg ti (n). Thus the motion of the quasi-particle representing an extra- 
charge in the disordered crystal is controlled by an Anderson hamiltonian with correlated 
diagonal and non-diagonal disorder that we shall now study in more detail. 



\ 
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FIG. 1: Distribution of molecular dipoles induced by a charge situated at the origin within the 
(a, b) plane. The arrows are projections of the dipole vectors on the plane. 

For this purpose we have developed a numerical method for calculating E p and J in 
disordered situations. Still working within the dipolar approximation, we solve the linear 
equations in relation [TJU numerically. This calculation is precise enough to determine the 
characteristics of the electronic polaron in naphtalene and anthracene. In tetracene and 
pentacene where the molecular polarizability is higher, corrections beyond the dipolar ap- 
proximation are important. In this case, the linear calculation would give values of X^i(n) 
exceeding 1/2 on certain sites, which would correspond to values of the dipole components 
exceeding the corresponding dipolar transition component /i^ of relation ^ In order to 
limit this artefact and to take into account the hyperpolarizability, the terms e |X^j(n)| 2 in 

equation El are replaced by - ^1 — y/l — A\X^i(n)\ 2 ^ . The system of equations in ITU1 then 

becomes nonlinear and must be solved by using an iterative numerical procedure. Figure Q 
represents the dipole distribution induced by a charge in a perfect crystal of pentacene. 

We have also introduced an orientational disorder which is known to be important in 
acenes. Pope and Swenberg give typical thermal disorder values for rotation angles of ±3 
degrees at room temperature [17]. Other sources of rotational disorder with angles even 
larger than 3° are also likely in actual acene materials, for example, metastable polymorph 
domains, stacking faults, zones close to grain boundaries, dislocation cores etc. |l8| . The 
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TABLE I: Polarization energy E p and renormalized transfer integral J/ J in the perfect crystals of 
naphtalene, anthracene, tetracene and pentacene for an angular disorder of 3°. 





Nph 


Ac 


Tc 


Pc 






21,5 


33,9 


48,2 


91 


Polarization tensor a 


(A3) 


17,6 


29,2 


34,7 


38 






10,1 


12,9 


15,6 


17,7 


Polarization energy! 20 ! E p 


(eV) 


-0,99 


-1,19 


-1,39 


-1,55 


Bare tranfer integral^ 21 ' J 


(meV) 


41,5 


47,9 


68,8 


97.8 


Renormalized transfer J/ J 


0,86 


0,83 


0,74 


0,64 


Diagonal disorder 5 


(meV) 


6,2 


9,7 


11,3 


15,8 


Correlation factor 5' 


(eV) 


5,3 


3,4 


1,0 


0,4 



simplest disorder distribution that is worth investigating is an angular gaussian distribution 
with mean square value A6. By averaging over 4000 samples for each A6, the mean square 
fluctuations 

5 = y/(E p 2 - <£p> 2 > 

and 8' can be deduced as functions of A6. Table U summarizes the results of our calcula- 
tions. Polarization energy E p and renormalized transfer integral J/ J in the perfect crystal 
are calculated from equations [TT] et [13] respectively. The diagonal energy disorder S is cal- 
culated for an angular disorder of 3° while 5' represents the correlations between diagonal 
and nondiagonal disorder (see equation 115)) . The polarization tensors and the bare transfer 
integrals are from litterature. (IlJ 

On the basis of this result, a good approximation for the related Anderson hamiltonian, 
would be 

H = J2 5 n\ n ) (n\ ~ z {6n+5n+h)l& ' \n) (n + h\ (15) 

n n,n+h 

where the distribution of 5 n is characterized by the mean square value 5 and the corre- 
lations by the value 5'. In the two dimensional channel of an organic field effect transistor 
these disorder energies, reported in table |H are non-negligible fractions of the renormalized 
transfer integral J . 

It is important to note that from the point of view of the extra-charge traveling into 
the molecular semiconductor, any thermal or static disorder which modulates the bare tight 
binding parameters is amplified by molecular polarization. These solvation effects yield much 
larger disorder energies than the corresponding single-electron quantities. The amplification 
factor is of the order E p /J ~ 10 (claim (ii) above). Moreover the electronic polaron effect 
introduces long-range correlations which are very efficient in terms of localization (claim 
(Hi) above). Work on the spectrum of hamiltonian is presently in progress to examine 
in detail the consequences of these correlations. At the present stage, we believe that the 
possibility to map the many-body hamiltonian El on the single electron hamiltonian El m 
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the subspace of the class of pertinent wave functions, will be useful for further theoretical 
developments in the field. 
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